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I. INTRODUCTION 



The probability density function (PDF) of velocity fluctuations is studied experimentally for grid 
turbulence in a systematical manner. At small distances from the grid, where the turbulence is 
still developing, the PDF is sub-Gaussian. At intermediate distances, where the turbulence is fully 
developed, the PDF is Gaussian. At large distances, where the turbulence has decayed, the PDF is 
hyper-Gaussian. The Fourier transforms of the velocity fluctuations always have Gaussian PDFs. 
At intermediate distances from the grid, the Fourier transforms are statistically independent of each 
other. This is the necessary and sufficient condition for Gaussianity of the velocity fluctuations. At 
' £H ', small and large distances, the Fourier transforms are dependent. 
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While velocity differences in turbulence u(x + Sx) — u(x) have attracted much interest, velocity fluctuations u(x) 
' themselves are also fundamental to describing the turbulence. Usually it is observed that the probability density 
function (PDF) of the velocity fluctuations is close to Gaussian g|. However, some experiments yield a sub- 
Gaussian PDF |3|, fj, pj, which has a less pronounced tail than a Gaussian PDF. The reason remains controversial 
because there have been no systematical studies. 

The observed Gaussianity had been explained by applying the cental limit theorem to the Fourier transformation of 
the velocity fluctuations jl] . This theorem ensures that a sum of many independent random variables has a Gaussian 
PDF However, owing to the steep energy spectrum of turbulence, the Fourier transforms have considerably 
different magnitudes. In this case, the central limit theorem is not applicable ||. We first reconsider the condition 
, for the velocity fluctuations to have a Gaussian PDF (Sec. II). 

■ For studies of velocity fluctuations, an experimental approach is preferable to the popular direct numerical simula- 
tion. Although the direct numerical simulation is useful in studying small-scale motions of turbulence, they are not of 
our interest because the velocity fluctuations are dominated by energy-containing large-scale motions. Those in the 
O [ direct numerical simulation suffer from its artificial boundary condition, initial condition, and forcing. 
^ ■ We obtain experimental data of velocity fluctuations in a wind tunnel downstream of a turbulence-generating grid 



(Sec. [II). The grid turbulence is not homogeneous in the mean- wind direction and thus allows us to study developing, 
tIh i fully developed, and decayed states of turbulence by increasing the distance between the anemometer and the grid. 
Oh- We find that the PDF of velocity fluctuations changes accordingly from a sub-Gaussian to Gaussian, and to hyper- 



Gaussian (Sec. IV). Here a hyper-Gaussian PDF has a more pronounced tail than a Gaussian PDF. We discuss the 



observed behavior by using the velocity fluctuations themselves, the velocity differences, and the Fourier transforms. 
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II. THEORY FOR GAUSSIANITY 



Suppose that velocity fluctuations u(x) are measured in a turbulent flow repeatedly over the range < x < L. The 
length L is much greater than the correlation length l u = J (u(x + 8x)u{x))d8x / (u 2 ) . Here (■) denotes an ensemble 
average. The measurements serve as realizations of the turbulence. Each of them is expanded into a Fourier series as 

u(x) = 22 a n{u) cos ( - J + b n (u) sin ( — J . (1) 

The Fourier transforms a n (u) and b n {u) for n ^> 1 have Gaussian PDFs over the realizations (see also Ref. |§| for 
a more mathematical explanation). This is because, for example, the Fourier transform a n (u) is obtained as 

»n(«) = y I u(x)cos( — - — ] da; 




...dx + ...+ / ...da: , (2) 

J(m-l)L/m J 

with 1 ^ m < n. The segment size L/m is set to be large enough so that the correlation (u(x + Sx)u(x)) has converged 
to zero at 5x — L/m. The integrations J^ m ... dx, •" ^ Xi > •••> an< ^ /(m-i)i/m "■ ^ are re g ar ded as independent 

random variables of the same magnitude. Then a n (u) has a Gaussian PDF as a consequence of the central limit 
theorem. The variance (a 2 (u)} is equal to the energy spectrum E n Q . 

The above discussion is not applicable to the Fourier transforms a n (u) and b n {u) for n ~ 1, i.e., those for large 
wavelengths. Nevertheless, as far as the wavelength is finite, it is possible to show Gaussianity of the corresponding 
Fourier transform by increasing the data length L and hence the n value. In the limit L — > oo, the transforms a n (u) 
and b n (u) for n ~ 1 become zero and do not contribute to the velocity fluctuations u(x). We are able to assume 
safely that all the Fourier transforms have Gaussian PDFs. Since this fact is independent of detailed dynamics, it is 
universal. 

If and only if all the Fourier transforms a n (u) and b n (u) are statistically independent of each other, the Gaussianity 
of them leads to Gaussianity of the velocity fluctuations u(x). To demonstrate this, we use the characteristic functions 
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Here E n cos 2 (2nnx / L) and E n sin 2 (2nnx / L) are the variances of a n (u) cos(2nnx / L) and b n (u) sin(2Trnx / L) . From 
the expansion formula (Q) and the independence of the Fourier transforms a n (u) and b n (u), it follows that the sum 

of the logarithms of the characteristic functions (f>^ (r) and <pn (t) is equal to the logarithm of the characteristic 
function 0(r) for the PDF of the velocity fluctuations u(x): 

oo 2 00 

ln0(r) = J2 M + = E E - ( 4 ) 

n— 1 ' n=l 

Thus the velocity fluctuations u(x) have a Gaussian PDF with the variance (u(x) 2 ) — X^Li ^n- The independence 
of the Fourier transforms also leads to the statistical independence of the velocity fluctuations u(x) at different spatial 
positions. 

Therefore, the necessary and sufficient condition for the velocity fluctuations to have a Gaussian PDF is the 
independence of the Fourier transforms. This is a good approximation for fully developed turbulence, where large- 
scale motions of energy-containing eddies are random and independent. Although turbulence contains small-scale 
coherent structures such as vortex tubes ||, their importance to the velocity fluctuations is negligible, i.e., being as 
small as the energy ratio of the dissipation range to the energy-containing range. 



III. EXPERIMENTS 



The experiments were done in two wind tunnels of Meteorological Research Institute. Their test sections were of 
0.8 x 0.8 x 3 and 3 x 2 x 18m in size (hereafter, respectively, the small and large tunnels). The small tunnel was 
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used to study developing and fully developed states of grid turbulence, while the large tunnel was used to study fully 
developed and decayed states. 

Turbulence was produced by placing a grid across the entrance to the test section. The grid consisted of two layers 
of uniformly spaced rods, the axes of which were perpendicular to each other. We used different grids in the small 
and large tunnels. The separation of the axes of their adjacent rods were 0.10 and 0.40 m, respectively. The cross 
sections of the rods were 0.02 x 0.02 and 0.06 x 0.06 m, respectively. The mean wind was set to be U ~ 10 ms -1 in 
the small tunnel and U ~ 20 ms -1 in the large tunnel. 

We simultaneously measured the streamwise (U + u) and transverse (v) velocities. They are velocity components 
that are parallel and perpendicular to the mean-wind direction, respectively. The measurements in the small tunnel 
were done on the tunnel axis from d — 0.25 to 2.00 m downstream of the grid with an interval of 0.25 m. Those in the 
large tunnel were done from d = 3.00 to 17.00 m with an interval of 1.00 m. The ranges of the measurement positions 
were to the limit of mechanical constraints of the wind tunnels. Since there was no overlap in the distance d between 
the small- and large-tunnel measurements, the individual data are identified by their d values. 

We used a hot-wire anemometer, which was composed of a crossed-wire probe and a constant temperature system. 
The wires were 5/xm in diameter, 1.25 mm in effective length, 1.25 mm in separation, and oriented at ±45° to the 
mean-wind direction. The wire temperature was 280°C, while the air temperature was 29-30°C in the small tunnel 
and 14-19°C in the large tunnel. We calibrated the anemometer before and after the measurements. 

The signal was low-pass filtered with 24dB/octave and sampled digitally with 16-bit resolution. In the small- wind 
measurements, the filtering was at 8 kHz and the sampling was at 16 kHz. In the large- wind measurements, the 
filtering was at 20 kHz and the sampling was at 40 kHz. The entire length of the signal was as long as 5 x 10 6 points. 
We obtained longer data of 2 x 10 7 points at the positions d — 0.25, 2.00, 8.00, and 12.00 m. 

The turbulence levels, i.e., the ratios of the root-mean-square values of the velocity fluctuations (u 2 ) 1 / 2 and (v 2 ) 1 ^ 2 
to the mean streamwise velocity U , were always low (< 0.2; see below). This is a good characteristic of grid turbulence 
and allows us to rely on the frozen-eddy hypothesis of Taylor, d/dt = — Ud/dx, which converts temporal variations 
into spatial variations in the mean-wind direction. 

The resultant spatial resolution is comparable to the probe size, ~ 1 mm. Since the probe is larger than the 
Kolmogorov length, 0.1-0.2 mm, the smallest-scale motions of the flow were filtered out. The present resolution is 
nevertheless typical of hot-wire anemometry [Q, f|, [l(J [fl|, [l2| . 

Fig. |Tj shows the mean streamwise velocity U, the root-mean-square fluctuations (u 2 ) 1 / 2 and (w 2 ) 1 / 2 , the correlation 
lengths l u and l v , the Taylor microscale A = ({u 2 } j '{(du / 'dx) 2 )) 1 / 2 , and the turbulence levels (u 2 ) 1 / 2 /U and {v 2 ) 1 / 2 j\J . 
For these and the other similar diagrams, open symbols denote the streamwise velocity while filled symbols denote 
the transverse velocity. The flow parameters change systematically with the distance from the grid, indicating a 
systematical change of the turbulence. This is especially the case in the turbulence levels. Our following studies suggest 
{u 2 } 1 / 2 /!/ > 0.1 for the developing state, (u 2 ) 1 / 2 /U ~ 0.04-0.1 for the fully developed state, and {u 2 ) 1 ' 2 ^ < 0.04 
for the decayed state. 

IV. RESULTS AND DISCUSSION 
A. Overview 

Fig. U demonstrates that Fourier transforms of the velocity fluctuations have Gaussian PDFs at the positions 
d = 0.25, 2.00, 8.00, and 12.00 m. The individual data were divided into 4864 segments of 2 12 points, which were 
regarded as independent realizations of the turbulence. They were windowed by a flat-topped function, which rises 
from zero to unity in the first small fraction of the data and falls back to zero in the last small fraction. The PDFs 
shown in Fig. |^ are those for the wave number k = n/L = 7 m -1 , where the energy spectra at d = 0.25 m have a peak 
(see below) . Since the PDFs of a n and b n should be the same, they were put together in order to minimize statistical 
uncertainties. We also obtained Gaussian PDFs of the Fourier transforms at the other wave numbers and at the other 
positions in the wind tunnels. 

However, velocity fluctuations do not necessarily have Gaussian PDFs. Fig. || shows the PDFs at the positions 
d = 0.25, 2.00, 8.00, and 12.00 m. The transverse- velocity PDF is sub-Gaussian at d = 0.25 m, Gaussian at d = 2.00 
and 8.00 m, and hyper- Gaussian at d = 12.00 m. The streamwise- velocity PDF tends to be skewed owing to a shear 
flow. 

Fig. |J shows the flatness factors F u = (u 4 ) / (u 2 ) 2 and F v , the skewness factor S u = (it 3 ) / (u 2 ) 3 / 2 , and the streamwise- 
transverse correlation C uv — ((u 2 v 2 ) — (u 2 ){v 2 ))/[({u 4 ) — (u 2 ) 2 )((v 4 ) — (w 2 ) 2 )] 1 / 2 . With an increase of the distance 
from the grid, the transverse- velocity PDF changes from sub-Gaussian (F < 3) to Gaussian (F = 3), and to hyper- 
Gaussian (F > 3). The streamwise- velocity PDF tends to be skewed (S ^ 0) and changes from hyper-Gaussian to 
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Gaussian, and to hyper-Gaussian. Also at large distances from the grid, the streamwise and transverse velocities have 
a significant correlation (C uv 3> 0). 

Since the streamwise fluctuations suffer from a shear, we are interested mainly in the transverse fluctuations. Their 
sub-Gaussian, Gaussian, and hyper-Gaussian PDFs are studied separately in the following subsections. 

B. Sub-Gaussian PDF in developing turbulence 

The transverse velocity has a sub-Gaussian PDF at the smallest distance from the grid, d = 0.25 m (Figs. || and 
^). Fig. |5|(a) shows energy spectra of the streamwise and transverse velocities. They have peaks at the wave numbers 
k ~ 7 and 14 m -1 . The flow is in a transition state from quasi-periodic motions due to wavy wakes of the grid rods to 
weak turbulence. Fig. Mb) shows the correlation coefficient between the Fourier transforms at adjacent wave numbers, 
Cnn' — ( a n a n')/(( a n)("n')) 1 ^ 2 wrt h n' = n + 1. At around the energy peaks, we observe significant correlations. The 
Fourier transforms are not mutually independent. 

The quasi-periodic motions should have finite correlation lengths. Namely, the motions should not be exactly 
periodic. If this were not the case, the central limit theorem would not be applicable to the calculation of the Fourier 
transforms (Sec. II). Then the transforms would not have a Gaussian PDF. This is inconsistent with our result 
(Fig-!). 

For velocity differences Su — u(x + Sx) — u(x) and Sv — v(x + 5x) — v(x), where Sx is the separation in the mean- 
wind direction, Fig. || shows the variances (Su 2 ) and {Sv 2 }, the flatness factors Fg u and Fg v , the skewness factors Sg u 
and Sg v , and the streamwise-transverse correlation Cg u g v . At small separations, we observe enhancements of Fg u , 
Fg v , Ssu and CguSv They are due to small-scale coherent structures ]|, ^|, |l2[ [l3| and are not of our interest. At 
large separations, we observe oscillations of (Su 2 ), (Sv 2 ), F$ u , F$ v , Sg u , Sg v , and CguSv Their wavelengths roughly 
correspond to the wave number k ~ 7 m" 1 of the spectral energy peak. The oscillations of (Su 2 ) and (Sv 2 ) are in 
phase, while those of Sg u and Sg u are ±90° out of phase and those of Fg u , Fg v , and Cs u Sv are 180° out of phase. 

Our results are explained if the velocity field is a superposition of a few quasi-periodic motions and a random 
background. It is actually possible to reproduce the oscillations and their phases qualitatively with a few sinusoidal 
functions satisfying the solenoidal condition and a random Gaussian noise. 

The quasi-periodic motions are predominant at local maxima of the variances (Su 2 ) and (Sv 2 ). Since the flatness 
factor Fsv is locally minimal and less than 3 there, the transverse- velocity amplitudes of those quasi-periodic motions 
lie in a limited range (see also Ref. This discussion would apply to the streamwise velocity as well. The 

oscillations of Sg u and Sg v imply the presence of quasi-periodic motions with different wave numbers that are coupled 
with each other (see also Fig. |(b)). 

The background flow is predominant at local minima of the variances. Since the flatness factor Fg u is greater 
than 3 there, the streamwise background tends to be intermittent. This tendency is not significant in the transverse 
background, which exhibits Fg v ~ 3 even at its local maxima. The slight enhancement of Cg u g v implies that the 
streamwise and transverse components of the background flow tends to have a correlation. 

Overall, the observed sub-Gaussianity of the transverse fluctuations is due to quasi-periodic motions that have 
finite amplitudes. The hyper-Gaussianity of the streamwise fluctuations is due to a background flow that tends to be 
intermittent. 



C. Gaussian PDF in fully developed turbulence 

The transverse velocity has a Gaussian PDF at intermediate distances from the grid, i.e., d > 1 m in the small 
tunnel and d < 10m in the large tunnel (Figs. 3 and 4). This is because turbulence is fully developed there. Wavy 
wakes of the grid rods have evolved to energy-containing eddies that are random and independent. The corresponding 
Fourier transforms are thus independent. Actually we did not find spectral energy peaks such as those observed at 
d = 0.25 m. The wave-number-wave-number correlation in the energy-containing range k < 10m _1 is absent within 
the statistical uncertainty, \C nn < \ < 0.1. 

The grid turbulence at the intermediate distances does not exhibit large-separation oscillations of (Su 2 ), (Sv 2 ), F$ u , 
Fg v , Sgu, Sg v , and Cg u g v . The flatness factors Fg u and Fg v as well as the skewness factors Sg u and Sg v at large 
separations are identical to the Gaussian values. The streamwise-transverse correlation Cg u g v at large separations is 
absent. These are observed usually in fully developed turbulence |J [l0| [l3|| . 
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D. Hyper-Gaussian PDF in decayed turbulence 

The streamwise and transverse velocities have hyper-Gaussian PDFs at largest distances from the grid, d > 10 m. 
The streamwise-transverse correlation C uv is also enhanced there (Figs. |]and ||). For the separations Sx = 0.20 m ~ l u 
and 5x = 0.40 m ~ 2l u , Fig. ^ shows the flatness factors Fs u and F$ v and the streamwise-transverse correlation CguSv 
They increase with the distance. Since the grid turbulence has decayed (Fig. |l|), there remain only strong energy- 
containing eddies fl, which intermittently pass the probe. The enhancement of the flatness factors Fs u and Fg v is 
due to enhanced values of 5u and 6v, which should be associated with the surviving strong eddies. The streamwise- 
transverse correlation CsuSv is enhanced if Su and Sv are enhanced simultaneously at the positions of those eddies. In 
this case, velocity fluctuations u{x) and v(x) are also enhanced at the eddy positions, resulting in the hyper-Gaussian 
PDFs as well as the significant streamwise-transverse correlation C uv . 

The presence of large-scale spatial structures implies a correlation among the corresponding Fourier transforms. 
Since the spatial structures consist of many Fourier modes, the correlation is not local in the wave number space. 
Although we failed to detect any significant correlation, we found moderate correlations with |C nn '| ~ 0.2 between 
many wave numbers in the energy-containing range k < 10m _1 . 

Roughly at the position where the flatness factors F u and F v begin to differ from the Gaussian value, the skewness 
factor for the streamwise velocity S u changes its sign (Fig. ||). This is probably because turbulence becomes weak. 
The positive skewness is more significant at a smaller distance from the grid and is attributable to the turbulence 
itself. On the other hand, the negative skewness is more significant at a larger distance. Using a low-pass filtering 
technique, we ascertained the presence of long- wavelength motions along the mean- flow direction (k <C 1 m _1 ). Their 
amplitude is larger and their PDF is more negatively skewed at a larger distance from the grid. These long-wavelength 
motions are attributable to effects of the wind tunnel, e.g., wall effects and deceleration of the mean flow. 

E. Small scale statistics 

Small-scale coherent structures such as vortex tubes in fully developed turbulence have attracted much interest 
H §, P, The statistics that are studied most often are the flatness factors Fg u /g x and Fg v /g x and the skewness 
factor Sg u /g x of the velocity derivatives du/dx and dv/dx. An increase of the microscale Reynolds number Re> = 
(it 2 ) 1 / 2 A/ ' v is known to cause the increase of Fq u /q x and Fg v /g x and the decrease of Sq u /q x ||. Here v is the kinematic 
viscosity. We briefly summarize overall trends of Re>, Fg u /g x , Fg v /g x , and Sg u /g x over distances from the grid. They 
are not of our interest but are expected to be helpful in a future experimental research. 

The microscale Reynolds number Re a decreases with the distance as shown in Fig. ||(a). The flatness factors Fg u /g x 
and Fg v /g x increase as shown in Fig. ||(b). The skewness factor Sg u /g x decreases as shown in Fig. ||(c). Thus the 
dependences of Fg u /g x , Fg v /g x , and Sg u /g x on Re a are apparently opposite to those mentioned above. This is because 
the turbulence state changes from developing to fully developed, and to decayed with an increase of the distance. In 
order to compare with other experimental data, the measurement is required to be done at a position where the grid 
turbulence is in the fully developed state. 

V. CONCLUSION 

The PDF of velocity fluctuations was studied systematically for grid turbulence. At small distances from the grid, 
where the turbulence is still developing, there are quasi-periodic motions having finite amplitudes, and hence the 
PDF is sub-Gaussian. At intermediate distances from the grid, where the turbulence is fully developed, motions 
of energy-containing eddies are random, and hence the PDF is Gaussian. At large distances from the grid, where 
the turbulence has decayed, there remain only strong eddies, and hence the PDF is hyper-Gaussian. The Fourier 
transforms of the velocity fluctuations always have Gaussian PDFs, in accordance with the central limit theorem. 
At intermediate distances from the grid, the Fourier transforms are statistically independent of each other. This is 
the necessary and sufficient condition for Gaussianity of the velocity fluctuations. At small and large distances, the 
Fourier transforms are dependent. 

Our result serves as an example that the velocity fluctuations could have a sub-Gaussian PDF if there exist strong 
finite-amplitude motions. We suspect that this is the case in experiments of Noullez et al. || and Sreenivasan and 
Dhruva [{§. They obtained F v a 2.85 in a free air jet and F u = 2.66 in an atmospheric boundary layer, respectively 
(see also Refs. O, ^2|), for which no explanation has been proposed yet. The signal could suffer from finite-amplitude 
motions generated by the jet nozzle or the surface. It is of interest to analyze such data in the same manner as in our 
present work. 
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FIG. 1: (a) Mean streamwise velocity U and root-mean-square values of velocity fluctuations (u 2 ) 1 ^ 2 and (v 2 ) 1 ^ 2 . (b) 
Correlation lengths l u and l v , and Taylor microscale A. (c) Turbulence levels (u 2 ) 1 ^ 2 jU and (v 2 ) 1 ^ 2 /U . The abscissa is the 
distance d from the grid. The open circles denote the streamwise velocity u, while the filled circles denote the transverse 
velocity v. The horizontal dotted lines separate the turbulence levels for which the turbulence is developing, fully developed, 
and decayed. In the calculation of the Taylor microscale, the velocity derivatives were estimated as, e.g., du/dx = [8m(x + A) — 
9,u(x — A) — u(x + 2A) + u(x — 2A)]/12A, where A is the sampling interval. 



FIG. 2: PDFs of Fourier transforms of velocity fluctuations at the positions d — 0.25, 2.00, 8.00, and 12.00 m. The wave 
number k is 7m~ 1 . We vertically shift the PDFs by a factor of 10 3 . The open circles denote the streamwise velocity u, while 
the filled circles denote the transverse velocity v. The solid lines denote Gaussian PDFs with zero mean and unity standard 
deviation. 

FIG. 3: PDFs of velocity fluctuations at the positions d = 0.25, 2.00, 8.00, and 12.00 m. We vertically shift the PDFs by a 
factor of 10 3 . The open circles denote the streamwise velocity u, while the filled circles denote the transverse velocity v. The 
solid lines denote Gaussian PDFs with zero mean and unity standard deviation. 

FIG. 4: (a) Flatness factors F u and F v of velocity fluctuations, (b) Skewness factor S u . (c) Streamwise-transverse correlation 
Cuv The abscissa is the distance d from the grid. The open circles denote the streamwise velocity u, while the filled circles 
denote the transverse velocity v. The horizontal dotted lines indicate the values expected for independent Gaussian PDFs. The 
skewness factor for the transverse velocity S v is close to zero within the statistical error of about ±0.02. 

FIG. 5: (a) Energy spectra of velocity fluctuations at the position d — 0.25 m. (b) Correlation between Fourier transforms at 
adjacent wave numbers C nn i with n' = n + 1. The abscissa is the wave number k = n/L. 

FIG. 6: (a) Variances (Su 2 ) and (5v 2 ) of velocity differences 8u and Sv at the position d = 0.25 m. (b) Flatness factors F$ u 
and Fsv (c) Skewness factors Ss u and Sg v . (d) Streamwise-transverse correlation Cs u Sv The abscissa is the separation Sx. 
The horizontal dotted lines indicate the values expected for independent Gaussian PDFs. 
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FIG. 7: (a) Flatness factors F$ u and Fg v of velocity differences 5u and <5v for separations Sx = 0.20 m (circles) and 0.40 m 
(triangles), (b) Streamwise-transverse correlation CsuSv The abscissa is the distance d from the grid. The open symbols denote 
the streamwise velocity u, while the filled symbols denote the transverse velocity v. The horizontal dotted lines indicate the 
values expected for independent Gaussian PDFs. 

FIG. 8: (a) Reynolds number ReA. (b) Flatness factors Fg u /g x and Fg v /g x of velocity derivatives du/dx and dv/dx. (c) 
Skewness factor Sg u /g x . The abscissa is the distance d from the grid. The open circles denote the streamwise velocity u, while 
the filled circles denote the transverse velocity v. The skewness factor for the transverse velocity Sg v /g x is close to zero within 
the statistical error of about ±0.0f . 
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